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Abstract 



We study the spectrum and magnetic properties of double quantum dots 
in the lowest Landau level for different values of the hopping and Zeeman 
parameters by means of exact diagonalization techniques in systems of 
N=6 and N=7 electrons and filling factor close to 2. We compare our 
results with those obtained in double quantum layers and single quantum 
dots. The Kohn theorem is also discussed. 
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I. Introduction 



The understanding of the structure and properties of double quantum dots (DQD's) 
grown on the perpendicular direction to their plain, and subject to an external constant 
magnetic field B, has attracted special interest due to the wealth of new quantum states 
that are actually realized. This wealth does not only refer to the ground states (GS's), 
which provide a quite intricate phase diagram^ but also to the variety of different 
excited states. This allows to have a large quantity of different sets of properties of 
the DQD system that can model for instance, point contacts within a device built for 
electron transport. 

Special interest has the regime for which the magnetic field is so strong that only 
the lowest Landau level (LLL) is occupied but not strong enough that prevent any spin 
polarization. We will always assume that the second Landau level is far enough that we 
can ignore any mixture between Landau levels. In the symmetric gauge, the projection 
of the DQD Hamiltonian to the LLL is given by^ 

H = aM + (3N- A,S, - AtX + H^^t (1) 

where 

(2) 
and 

P-^yf4 + ^o, (3) 

ujc = is the cyclotron frequency, m* being the effective electron mass in the semi- 
conductor host and e and c are the electron charge and the speed of light in vacuum 
respectively. The frequency associated with the parabolic confining potential in both 
dots is given by a;o, M is the total angular momentum and N is the total number of 
electrons. The Zeeman coupling is given by = qhbB with g being the Lande ^^-factor 
and being the Bohr magneton ( for the free electron mass, Hb = eh/2mc). The single 
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Uj'i + Aujn — UJr 
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particle energy gap between symmetric (s) and antisymmetric (a) states, combinations 
of |r) and |/) single particle states of electrons confined in the right and left dot respec- 
tively (see below) is given by 2At and X — Ng — Na is the balance between symmetric 
and antisymmetric single particle states. Finally, Hint is the Coulomb interaction term. 
Hereafter all distances will be given in units of the magnetic length defined as: 

/ [ 

^ [m*(a;2 + 4^2)1/2 

and all energies in units of /{eIb) being e the dielectric constant of the semiconductor 
host. The Zccman energy, the tunneling term scaled by A^, the kinetic contribution 
given by a and /3 and the Coulomb interaction, provide the ingredients of the system. 
The applied magnetic field can be directed in any direction in space in such a way that 
its action has different effect in the kinetic contribution in which only the component 
along the Z-direction plays a role and in the Zeeman effect in which the total magnetic 
field contributes. The eigenstates of the Hamiltonian are characterized by the total 

— * 

angular momentum M, the total spin along the S-direction given by Sz and the parity 
P^'^ . These parameters are related with the invariance of the Hamiltonian under space 
rotations along the Z-direction, rotations of spin and specular reflection in the plane 
inter-dots respectively^. We will denote by (M, Sz, P) the configuration that determines 
each separated subspace of eigenstates. 

Previous studies in DQD's based on exact diagonalizations on one hand^'^ and the 
large experience extracted from mean field (MF) approximations and effective field the- 
ories in double layer (DL) systems on the other hand^~^, have provided a quite complete 
picture of the GS phase diagram. It is well estabhshed in both DL (Ref. 7) and DQD 
(Ref. 1) systems that the variation of , leaving fixed all the other parameters, induces 
several changes in the GS. For A^ ^ A^ the GS is ferromagnetic |FM) with symmetric 
and antisymmetric single particle states equally populated and with all the electrons 
in spin up states. In contrast for A^ ^ At the GS is symmetric |SYM) with all the 
electrons occupying symmetric single particle states with up and down spins equally 
populated. In between, for comparable values of Zeeman and tunneling energies, inter- 
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(4) 



mediate states of the type called canted states (|C)) become GS's. The canted state is 
the GS solution between the |FM) and |SYM) in a MF calculation for a DL system and 
it is characterized by a ferromagnetic order in the direction perpendicular to the layers 
and antiferromagnetic order in plane^. 

A remarkable result of the Hartree Fock (HF) approximation for DL quantum Hall 
systems at fiUing factor u — 2 is that there are only three possible different GS's: |FM), 
|G) and |SYM) distributed over a universal phase diagram which boundaries depend only 

on three energy scales whereas, the Hamiltonian depends on four independent energy 
scales^. 

Certain ground states obtained from exact diagonalization in DQD's have also been 
identified as canted states since their overlap with the MF canted GS's for finite systems, 
projected on the appropriate subspace with well defined S^, is nearly one^. The total 
number of different canted states between |FM) and |SYM) in a DQD depends on the 
number of electrons. As a general rule, keeping fixed and increasing A^, the GS 
evolves as: 

I FM) (M, S,, P) (M, - 1, -P) (M, 5, - 2, P), ... , I SYM) (5) 

A simultaneous change of spin and parity occurs for the GS in each phase transition 
whereas the total angular momentum remains unchanged^. 

No direct information about the order in plane has been obtained so far from the 
exact solutions due to the fact that they have all the quantum numbers (M, S'z, P) well 
defined. No information can directly be obtained from order parameters that depend on 
single operators of the type or since their expectation values vanish. However, 
indirect information about the order in plane can be obtained from the properties of 
the excited states as we will show in the discussion of Figs.(l) and (2) below. Within 
the great number of possible low energy excitations that provide information about the 
properties of the system, we will draw our attention to those related to three main 
different points: 
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First, we search for particle-like excitations independent of the electron-electron in- 
teraction and so independent of the number of electrons. We discuss when the Kohn 
theorem^° holds in a parabolic DQD system ( as it is always the case for a single parabolic 



Second, we look for the types of excitations that soften as they come close to a GS 
transition that takes place with the variation of some input parameter. These excitations 
provide a clear and easy way to map the GS phase diagram. 

Third, in the last point we concentrate in the dispersion relation u!{l) of two different 
type of excitations, one with and the other without a simultaneous spin and isospin 
flip (see below). In a single QD the Coulomb contribution to the dispersion relation of 
excitations over a ferromagnetic GS decreases with angular momentum due to expansion. 
This is the general behavior, except at the "magic" values of l*^^ for which the system 
increases in angular momentum from 1* to T -|- 1 but does not increase in size and hence 
the Coulomb contribution remains constant. These magic values of I are related to the 
incompressible states. Our aim is to see if there is a similar behavior in the DQD system. 

This paper is organized as follows: In Section II we make a brief account of the 
exact diagonalization used in our calculations. In Section III we present the outstanding 
features of some excited states of different multipolarity over different GS's, paying 
special attention to the different behavior of even and odd systems. We will follow the 
three main points previously mentioned. Finally in Section IV we draw our conclusions. 



Exact diagonalization can be performed in separate boxes for different conflgurations 
characterized by (M, Sz,P). Each configuration determines a finite-eigenstate subspace. 
Within the LLL regime, we will use in all expansions the non-interacting single particle 
wave functions which do not have nodes along the radius and are given by 



QD) 



II. Exact diagonalization 




(6) 
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where m is the single particle angular momentum, A — r,l, cr — ^, i, f—{x,y), and 
(prni'f^ are the Fock Darwin wave functions given by^^ 

^ V27rm!2'" ^ ^ 

Along the Z-direction we assume delta charge distributions separated by d, the distance 
between the dots (we consider o? = 1 in all numerical performances). 

Within our calculations we will use symmetric |s) and antisymmetric \a) single par- 
ticle states, related by |r) and |/) by 

I s) = -^i I r) + I ) (8) 



and 
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r)- 10 ), (9) 



we will use the term isospin referred to this degree of freedom (sometimes referred as 
pseudospin in literature, see Ref.^). None of the parameters r and I or s and a are well 
defined quantum numbers since the Coulomb interaction mixes s and a and the tunneling 
process mixes r and /. The well established restriction is that parity must be preserved, 
this means that the change in symmetry ( s — > a or a — > s) due to electron interaction 
must take place by pairs of electrons, and never by one alone (see below). 

We will consider excited states over the three possible types of GS's and choose the 
parameters in such a way that for an even number of electrons in the ferromagnetic 
GS we have filling factor u — 2. We proceed as follows. Once the input parameters 
are fixed, we determine first the finite number of Slater determinants built up from 
single particle wave functions of the type given by Eq.(6) which provide a bases for 
each subspace configuration {M,Sz,P). Then, the diagonahzation of the Hamiltonian 
given by Eq.(l) is straightforward except for the Coulomb term. Although the Coulomb 
interaction docs not mix |r) and |/) single particle states, some manipulations must be 
done in order to operate over \s) and \a) wave functions. In second quantization, the 
interaction Hamiltonian is given by 



Hint = ^ Yli^i^j I ^ I (t>k(t>i)(4a'^aiak (10) 

ijkl 

where V is the Coulomb interaction and the sub-indexes denote angular momentum 
m — 0,1,2,.., spin= |, J, and isospin t— a, s. Taking into account that the single 
particle wave functions are related by Eq.(8) and (9) and that the right and left single 
states are ^-distributions of the type 5{z) and 5{z — d) respectively, it is easy to show 
that 



^lct>a = <P:<Ps 1{S{Z) +d{z-d)) (11) 



and 
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~ -{8{z) - 5{z - d)) (12) 



for each electron. As a consequence, there are only three possibilities for the expectation 
values of V: (i) the interacting electrons do not change their isospin as in 

(0,0, I V I (t>s(t>s) = \{ {Vrr) + {Vrl)) (13) 

(ii) only one electron changes its isospin as in, 

{ct>sCt>a\V\<i>aCt>a)=Q (14) 

and (iii) both electrons change their isospin index as in 

{<t>s<i>a I V I = ^{{Vrr) - (F^)) . (15) 

In the brackets on the right hand side of these equations, the integral over the Z- 
coordinate has been performed and the potentials are given by 

2 2 
Vrr ^Vu^-. Vrl = V,r = , , ^ (16) 

with r —\ r*! — r2 |, 1 and 2 being the two interacting electrons. That is to say, the 
Coulomb interaction either leaves unchanged the isospin of the electrons, if it operates 
with 
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Vo = \{Vrr + Vrl) (17) 

or changes the isospin of two electrons, if it operates with 

= ^(K. - Vrl) (18) 

in such a way that parity P, given by P = (—1)"^/^ is preserved {X = Ng — NaY ■ The 
change of isospin of a single electron is forbidden. 

Ill Multipolar excitations 

Guided by an interest in the properties of the multipole mode excitations and also 
by an interest in those excitations that soften at the boundaries of the phase transitions 
within the phase diagram, we explored the eigenenergies and eigenstates coming out of 
exact diagonalizations. We consider, in general, the multipolar excitations characterized 
by uj{l) with or without a simultaneous change of spin and/or parity and in some cases 
we concentrate on I — Q ot I — 1. 

Before focusing on the first point, that is, within the search of excitations which keep 
the internal Coulomb energy constant, let us briefly comment on the Kohn theorem for 
a DQD. The Kohn theorem for a single QD states^^ that in the process of absorption (or 
emission) of a photon of wavelength much larger than the radius of the dot confined by 
a parabolic potential, the initial and final electronic states can differ from one another 
only by the center of mass (CM) excitation. The number of confined electrons and the 
interaction between them has no influence on the values of resonance energies. For a 
DQD, however, qualifications to this statement may be required for different directions 
of the incident electric field with respect to the plane of the dots. It is easy to see that for 
homogeneous electric fields in-plane the Kohn theorem also holds for DQD's. However, 
this is not the case for homogeneous electric fields along the Z-direction, as we argue 
below. 

Taking into account Eqs. (8)- (9) and the expression of the exciting potential (asso- 
ciated with E±_ directed along the Z-direction) 
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HE^^Y^^i^^ ' (19) 



i=l 



it is easy to see that (0^ | He^ \ (pa) ^ d ^ 0. That is to say, the operator Hej^ changes 
the parity of the system and so [ Hej^ , Hint] 7^ since the eigenstates of Hint have well 
defined parity. The part of the Hamiltonian which depends on the relative coordinates 
only receives other contributions aside from the interaction Hamiltonian (coming from 
the kinetic and from the tunneling tcrms^^ within the first quantized expression) however 
the possibility that the last bracket can be compensated by the brackets of the other terms 
can be disregarded as they depend on different independent parameters. An alternative 
reasoning is presented in Appendix A. In the analysis of the exact diagonalization results 
which follows, we will verify, among other things, the applicability of the Kohn theorem. 

Let us begin with the N=7 electron system and consider first the excitations of 

the |FM) GS that increase the angular momentum in one unit and leave all the other 
parameters of the configuration unchanged, that is to say, we consider, 

{M,S,,P) {M+1,S,,P). (20) 

The system jumps from a 1-dim space configuration to a 2-dim space. The result is 
the excitation of the CM by a leaving the internal Coulomb energy unchanged. This 
is the well known intra-Landau level dipole excitation whose energy decreases when 
the magnetic field increases that is, the — a far-infrared resonance (FIR)^^. No 
inter-Landau level transition of energy given by 
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(21) 



is possible within the LLL considered in our calculation. Other Coulomb-independent 
interaction excitations are possible in the N=7 system. If we also allow changes in parity, 
that is to say, if we consider, 

(M, S,, P) ^ (M + 1, S,, -P) (22) 

the CM is excited with an energy a -|- (it is a CM excitation since the same result would 
be obtained for N=l). The previous case (of energy a) would correspond to the excitation 
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made by a nearly constant electric field ( we are assuming dipole approximation) directed 
along the X direction (£^11 ) and the last one (of energy a + would correspond to an 
electric field with an additional non- vanishing component along the Z direction (let us 
call it Exz)- Incidentally, this result can be used to determine experimentally the relative 
orientation of the DQD respect to the incident beam. The system absorbs a photon of 
energy a + At with maximum probability when the angle between the direction of the 
incident electric field and the normal to the plane of the DQD is ^ = 7r/4. 

For both excitations, with and without parity change, the system goes from a 1-dim 
space (|FM)) to a 2-dim space of excited states. The CM excitation leaves the system 
in the higher energy state within the excited 2-dim configuration in such a way that, 

E2{M + 1,S,,±P) = E,{M,S,,P) + a + (1^1)^2 (23) 

and 

Ei{M + l,S,,±P) < E2{M + 1,S,,±P) (24) 

where the cigcncncrgics Ei have been ordered from lower to higher energy within each 
subspace. See Table (cases A and B). The difference in energies in the excited config- 
urations come from the decrease of Coulomb interaction due to the expansion of the 
eigenstate of lower energy. The highest energy eigenstate is a compact state and does 
not change in size. The two excitations of case A with and without Coulomb contribu- 
tion are equivalent to the "sum mode" and the "difference mode" found by Girvin and 
MacDonald^ for DL (without tunneling or Zeeman terms and for filling factor u — 1/2) 
within a single-mode approximation. Once we know the decrease of energy of Coulomb 
origin from ^2(10, 7/2,0) - Ei(10, 7/2,0) = E^ = 0.0998, one can distinguish in the 
(11, 7/2, 0) configuration the double dipole Eq = Ei{9, 7/2, 0) + 2q; = 16.0802 (a = 0.2) 
from a quadrupole excitation given by E4 — Ei(9, 7/2, 0) + 2a — E^ — 15.9803. Surpris- 
ingly, the appearance of this cigenenergy in the output means that an eigenstate which 
is not a compact state nor a fully expanded state, is a possible realization. The full 
expansion would imply a larger value of Ec as AM = 2. 
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However, not all the previous scenario is reproduced in the N — 6 system. If by 
analogy we consider the excitations (M, S^, P) — * (M + 1, Sz, ±-P) from the |FM) GS, 
the energy E2{M + 1, Sz,P) — Ei{M, Sz, P) + ct is obtained according to the behavior 
of the N — 7 system but no trace of the Aj energy (At = 0.06) is present in E2{M + 
l,Sz,-P). 

Prom the analysis of the results of the N — 7 system, one would suspect that the 
CM excitations (of energies a or a + At) are a consequence of the Kohn theorem in the 
DQD system, in agreement with our assumption of a parabolic confining potential and 
dipole-excitation approximation. However, if it would be the case, the A'" = 6 system 
would behave in the same way. The reason for the difference is provided by the different 
action of the electric fields that excite the system. The field in plane produces a global 
shift of the system which does not affect the internal electron-electron interaction and 
the CM absorbs the total energy a. In contrast, if the external field has non- vanishing 
component along the Z-direction, aside of the global absorbtions of energy a an extra 
amount of energy is involved in a tunneling process. In the odd electron system, the 
unpaired electron tunnels from one dot to the other. This change, however, does not 
modify the electronic distribution asymmetry which was already present in the GS and 
consequently does not modify the Coulomb interaction. On the other hand, in the even 
electron system due to the electronic jump, a change is produced between a symmetric 
distribution to a non-symmetric one giving necessarily Coulomb contribution to the final 
energy. That is to say, the excitations induced by E± on an odd number of electrons 
does not modify the Coulomb energy. However this is not a consequence of the Kohn 
theorem, as one may erroneously conclude. 

Within the second main point, the search of excitations that mark the phase transi- 
tions, we concentrate on the evolution of the system as Aj increases, leaving all the other 
parameters unchanged. The appropriate excitations turned out to be the lowest energy 
excitations. Figs.(l) and (2) show the excitation energy as a function of Aj for the four 
different GS phases existing along the At energy variation line for iV = 6 and N = 7. 
In all the cases the angular momentum is preserved and a simultaneous spin and parity 
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flip takes place. That is to say, 

u = E^{M,S,±l,-P) - E,(M,S,,P) (25) 

where Ei{M, S^, P) is the GS energy at Aj and Ei{M, ± 1, —P) is the lowest eigenen- 
ergy within the excited conflguration at the same A^. These excitations provide the 
lowest value ( AM = / = ) of the multipole dispersion relation u}{l) given below. 

Figs. (3) and (4) show the dispersion relation of the two different types of multipole 
modes uj{l) mentioned previously within the third main point. In both cases, the lowest 
energy eigenstate within each subspace is considered. The values of Aj were chosen in 
such a way that for the two systems the GS is the first canted state: (6, 2, 1) for N=6 
and (9, 5/2, 1) for N=7. The general trend is given by the dashed fine which corresponds 
to, 

u = Ei{M + l,S„P) - Ei{M,S„P) (26) 

The increase of angular momentum in one unit increases the kinetic energy of the system 
by q; (a = 0.2). This increase is partially compensated by the decrease of Coulomb 
energy due to expansion. For both systems, the Coulomb contribution is numerically, 
for the particular values of the parameters that we have taken, about 2a/ 3 in such a 
way that the total amount of energy gained in each step is about a/3. This gives a quite 
monotonous behavior. 

In contrast, unexpected features were obtained for excitations of the type (solid line): 

oj = E,{M + l,S, + 1, -P) - E,{M, S,, P) (27) 

where spin and parity flip takes place simultaneously. For low values of /, the contribution 
of the Coulomb energy to uj is given by a more or less constant amount of, numerically, 
a/2. However, after several steps, 4 for N=6 and 6 for N=7 the system suffers a sudden 
expansion that reduces drastically the Coulomb interaction. In the N=7 case, due to 
the presence of an unpaired electron, the amount of the Coulomb contribution is so 
large that the total energy decreases. Looking at the occupancies of the single particle 
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states of the excited systems the nature of the sudden change becomes clear. For the 
ferromagnetic excited states of the type (6 + /, 3, 0) for N=6 and (9 + /, 7/2, 0) for 
N=7 only spin- up states can exist. The occupied single particle states are |st) and |at). 
As I increases from I — 0, there is a slow transfer of electrons from \a\) to |st)- The 
sudden change takes place when the following structure is possible: N-1 electrons in \s]) 
with m = 0, .., iV — 2 and one electron in \a'\) with m=0, which means M=10 (6+4) for 
N=6 and M=15 (9+6) for N=7. For this eigenstate almost all the weight is in a single 
Slater determinant which means that the system is nearly uncorrelated. 

The unexpected result is in clear contrast with the result obtained for a single QD at 
the magic values of the angular momentum, where all the energy is absorbed by the CM. 
In the DQD most of the absorbed energy is transformed into internal energy releasing 
the electrons from their interaction. 

IV Discussion and Conclusions 

We have analyzed several types of low energy excitations over the three possible GS 
of a DQD confined by parabolic potentials in each plane and separated by a distance d. 

The LLL regime was considered and the input parameters were chosen in such a way 
that the filling factor of the ground states and some of the excited states is close io u — 2. 

From the study of several dipole excitations ( AM = / = 1) we verify that the Kohn 
theorem is preserved in a DQD if the incident homogeneous electric field is directed along 
the plane of the dots. For systems with an odd number of electrons, additional excitations 
independent of the Coulomb interaction, as the one induced by are possible. This 
last possibility is due to the fact that in this case, the isospin change for a single electron 
does not change the Coulomb energy of the system, although it is not a consequence of 
the Kohn theorem. 

The lowest energy excitations arc characterized by a simultaneous spin and parity 
flip keeping the angular momentum unchanged ( spin-density-waves (SDW's) for / = 0, 
with energy defined by u see Figs. 1 and 2). They turn out to be the appropriate 
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excitations which become degenerated with the GS at the transition boundaries in the 
phase diagram. 

Due to the extra degree of freedom (as compared with a single dot or a single layer), 
represented by the isospin states and due to the interplay between tunneling and Coulomb 
interaction, the energy lu has a Coulomb contribution. This is in contrast to the case of 
a single layer for which, in the limit A; — > the non-interacting contribution given by 
Az is recovered. 

The softening of the cu modes with the variation of shows phase instabilities previ- 
ously detected in the determination of the GS phase diagram from exact diagonalization^. 
There are some similarities and some differences between our results and those reported 
by Das Sarma et al.^ within a HF calculation for a double layer system. We will fol- 
low the arguments given by these authors to analyze our results. Within the |FM) and 
|SYM) phases, the structure of the curves co/At is similar, in both cases the u; mode 
softens as it approaches the phase transition from |FM) to |C) and from |SYM) to |C) 
respectively. In addition, the increase of u as it moves away from the boundary is larger 
in the symmetric phase in both calculations. However, within the canted state we obtain 
finite values of u; although much smaller than those in the |FM) or |SYM) phases, in 
contrast to the results obtained in the DL for which a; = over the full canted phase. 
Aside from the previous outcomings, there is another main difference: the canted state 
in our calculation is an eigenstate of the Sz operator whereas this is not the case for 
the HF canted phase. At the boundaries separating different phases, however, due to 
energy degeneracy, the superposition of states of different well defined spin gives rise to 
states which are not eigenstates of S^. Das Sarma et al. have proved that the existence 
of a gapless mode is directly due to the canted antiferromagnetic spin ordering. As a 
consequence, even though the canted states in a DQD may be interpreted as having 
antiferromagnetic order in the plain of the dots^, the order is not complete, producing 
gapped excitations probably due to edge effects. The exceptions are at the boundaries 
separating different phases where due to degeneracy, the gapless mode is recovered. 
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Finally, unexpected nearly uncorrelated states have been found (see Figs. 3 and 4), 
produced by SDW excitations over the canted ground states, that develop a sudden 
expansion leaving the system in a ferromagnetic state close to 1/ — 1. Besides that, the 
structure developed by the sudden expansion on the dispersion relation at a particular 
value of I is understood as a consequence of the existence of an intrinsic length scale in 
the system, related to the antiferromagnetic spin ordering in plane. The structure is well 
defined in the SDW dispersion relation (solid line in Figs. 3 and 4) and only slightly 
apparent in the charge density wave curve (dashed line in Figs. 3 and 4) due to the fact 
that within the LLL, local variation of spin is coupled to variation of density^ producing 
a larger effect when spin and charge density fluctuations are induced simultaneously by 
the excitation. No structure was found in similar dispersion relations of SDW over the 
ferromagnetic or symmetric phases. In these last cases the dispersion is given by nearly 
parabolic smooth curves whereas, over the canted states, in the long- wavelength limit, 
the dispersion relation is found to be linear, consistent with the fact that it describes 
antiferromagnetic fluctuations. 

Within the HF calculation performed by Das Sarma for DL, the SDW dispersion 
relation does not have this type of structure. This difference may be related to the 
different role that the kinetic contribution to the energy plays in extended and in confined 
systems. 

Finally, in the hmit I — > 00 the dispersion relation of the SDW does not approach 
asymptotically a constant value as it is the case for a single layer (SL) or DL systems due 
to the fact that at this limit the excitation energy always recovers the single particle value. 
For SL or DL this energy is given by uu' + Vex where uu' is the non-interacting excitation 
value ( a combination of and A^) and Vex is the exchange single particle energy of 
an electron in the GS, independent of the linear momentum k due to the degeneracy 
existing in extended systems. In contrast, in a DQD the parabolic potential breaks the 
degeneracy producing an increase of energy with increasing angular momentum. This 
gives a nearly linear curve as Z — > 00 typical of a single particle. 
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No total spin or space correlations in the ground states have been investigated through 
the density-correlation functions which may signal additional symmetry breaking effects. 
Symmetry breaking effects of this kind have been found previously in literature for QD's, 
DQD's (vertical arrangement) and QD- molecules (horizontal arrangement) within un- 
restricted Hartree-Fock calculations for low magnetic fields^^'^^ and within Hamiltonian 
diagonalization for high magnetic fields^^. Whereas a comparison with the case of low 
magnetic fields is not possible because we have strong magnetic fields which project the 
system to the LLL, it would be interesting to see if our method, which includes spin de- 
grees of freedom and does not require any truncation, shows symmetry breaking effects 
similar to Ref.^^. This is left for future investigations. 

Acknowledgements 



We gratefully acknowledge C. Tejedor and L. Martin-Moreno for helpful discussions 
and for the code used for the Hamiltonian diagonalization. This work has been performed 
under Grants No. BFM2002-01868, No. FPA2001-3598 from MCyT and Feder (Spain) 
, and No. 2001SGR-0064 and No. 2001SGR-00065 from Generahtat de Catalunya . 



For simplicity we neglect spin indices since they are irrelevant for the following dis- 
cussion. We characterize the first quantized wave function of the N particle system as 
^ai,...,ajv (^1 ^n), where = 1,2 are layer indexes indicating the dot in which the 
i-th electron sits. Hence the first quantized Hamiltonian is a matrix in layer-index space 



An in-plane homogeneous electric field clearly produces a contribution proportional 
to the identity matrix in the layer space. 



and hence the interaction only depends on the center of mass coordinates, namely the 
Kohn theorem applies. A homogeneous electric field perpendicular to the plane E± 
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biai ;...;6iv^iV * 
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produces a contribution diagonal in layer space (but not proportional to the identity 
matrix) which reads 

i=l ^ \i=l 

where d is the distance between the two dots. 

The hopping term is not diagonal in layer space, 

N \ In 

Hhop ~ X! n ^bjaj Sb^ai JJ 
i=l \j=l J \j=i+l 

where Su — S22 — and S12 = S21 = 1- Then it is easy to check that 

^ i=l \3=l ) \j=i+l 

where e^^ai is the antisymmetric tensor. Hence a homogeneous electric field perpendic- 
ular to the plane changes the dynamics in layer space in a non-trivial way. Since the 
latter is entangled with the relative motion through the Coulomb term we conclude 
that a homogeneous electric field perpendicular to the plane may produce transitions 
between different states with the same center of mass quantum numbers, namely the 
Kohn theorem does not apply. 
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FIG. 1. Lowest energy excitation for A'^ = 6 as a function of for different regions of the 
phase diagram. The excited configurations are (from left to right) (6, 3, 0), (6, 2, 1) and (6, 1, 0) 
for the positive-gradient curves and (6,2,1), (6,1,0) and (6,0,1) for the negative-gradient 
curves. The values a = 0.2, = 0.02 and /? = 1.4 have been considered. 
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FIG. 2. The same as Fig.l for N = 7. The excited configurations are (from left to right) 
(9,7/2,0), (9,5/2,1) and (9,3/2,0) for the positive-gradient curves and (9,5/2,1), (9,3/2,0) 
and (9, 1/2, 1) for the negative-gradient curves. 
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FIG. 3. Dispersion relation io{l) of two different excited states: (6 + 1, 3, 0) (solid line) and 
(6 + 1, 2, 1) (dashed line) for = 6 and At = 0.07. The GS is the first canted state. The same 
values of a, P and as in Fig.l have been used. 
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FIG. 4. The same as in Fig. 3 for the excited states: (9+/, 7/2, 0) (solid line) and (9+/, 5/2, 1) 
(dashed line) for = 7 and A* = 0.085. 
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Table 1: Eigenenergies obtained from different configurations. We considered: 
At = 0.07, and At = 0.06, for N = 7 and = 6 respectively, in both cases 
a = 0.2, l3 = 1.4 and A^ = 0.02 and the GS is ferromagnetic. A is the excitation 
due to an electric field given hy E = E^^i, and B due to £" = E^y^ being i the 
unitary vector along the X-axis. C is due to an electric field expanded up to 
the quadrupolar approximation and directed along the X-axis. 



N=6 



GS 


A 


B 


G 


(6,3,0) 
El = 12.6802 


(7,3,0) 
El = 12.7868 
E'j = 12.8802 


(7,3,1) 
El = 12.7232 
E2 = 12.8434 


(8,3,0) 
El = 12.8804 
E2 = 12.8915 
E3 = 12.9438 
Ei = 12.9868 
^5 = 13.0802 


N=7 


GS 


A 


B 


G 


(9,7/2,0) 
El = 15.6802 


(10,7/2,0) 
El = 15.7803 
E2 = 15.8802 


(10,7/2,1) 
El = 15.8503 
E2 = 15.9502 


(11,7/2,0) 
El = 15.8700 
^2 = 15.8906 
Es = 15.9523 
^4 = 15.9803 
Es = 16.0231 
Ee = 16.0802 



23 



